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ABSTRACT 

We study the spectrum of high frequency radiation emerging from mildly dissipative photospheres of long- 
duration gamma-ray burst outflows. Building on the results of recent numerical investigations, we assume 
that electrons are heated impulsively to mildly relativistic energies by either shocks or magnetic dissipation 
at Thomson optical depths of several and subsequently cool by inverse Compton, scattering off the thermal 
photons of the photosphere. We show that even in the absence of magnetic field and non-thermal leptons, 
inverse Compton scattering produces power-law tails that extend from the peak of the thermal radiation, at 
several hundred keV, to several tens of MeV, and possibly up to GeV energies. The slope of the high-frequency 
power-law is predicted to vary substantially during a single burst, and the model can easily account for the di- 
versity of high-frequency spectra observed by BATSE. Our model works in baryonic as well as in magnetically 
dominated outflows, as long as the magnetic field component is not overwhelmingly dominant. 

Subject headings: gamma-ray: bursts — radiation mechanisms: non-thermal — methods: numerical — rela- 
tivity 



1. INTRODUCTION 

Understanding the origin of the prompt emission of 
Gamma-Ray Bursts (GRBs) has been hampered by two ma- 
jor challenges that have proven formidable obstacles despite 
more than four decades of studies. Radiating photons from 
an ultra-relativistic jet requires first a mechanism to convert at 
least some of the bulk kinetic energy into thermal motion of 
electrons, second a mechanism to radiate the electrons' en- 
ergy into photons. These two processes are known as the 
problem of the dissipation mechanism and of the radiation 
process. 

The standard model that has been developed over the years 
assumes that internal shocks are responsible for the dissipa- 
tion (Rees & Meszaros 1994), while electrons gyrating in a 
shock-generated magnetic field produce the radiation via the 
synchrotron process (Meszaros et al. 1994; Piran 1999; Lloyd 
& Petrosian 2000). Both mechanisms are fraught with numer- 
ous problems. 

Internal shocks assume that the outflow is released by the 
central engine with fluctuations in the Lorentz factor, so that 
different parts of the flow collide with each other, produc- 
ing shocks that dissipate energy. Unfortunately, the internal 
shock mechanism has very little predictive power, since any 
behavior of the light curve can be explained by a suitable 
choice of the ejection history of the central engine. Since 
the physics of the emission from the engine is largely un- 
constrained, it is very hard to disprove internal shocks based 
on any observation. The only robust prediction of internal 
shocks as a dissipation mechanism is the low efficiency of the 
process (Kobayashi et al. 1997; Lazzati et al. 1999; Spada 
et al. 2000), due to the fact that internal shocks can dissi- 
pate only the energy associated with differential motions and 
not the energy associated to the bulk motion, which is much 
larger. Observations, instead, show that at least for a fraction 
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of bursts, the efficiency of the prompt phase is 50 per cent, 
if not higher (Zhang et al. 2007). An alternative model that 
is becoming increasingly popular is magnetic dissipation in 
a Poynting dominated outflow (e.g., Thompson 1994; Spruit 
et al. 2001; Giannios & Spruit 2005). Even though magnetic 
dissipation could in principle provide very high efficiencies, it 
is, again, a very uncertain process and as such provides very 
little predictive power to allow for meaningful comparisons 
with observations. 

Synchrotron radiation, with possible self-Compton compo- 
nents, has long been considered the best candidate to explain 
the prompt emission of GRBs. While synchrotron radiation 
can easily explain the non-thermal nature of the observed 
spectrum, more thorough scrutiny reveals several important 
problems. First, optically thin synchrotron emission has a 
very well defined limiting slope a > -1 /3 (where F{v) oc j/"") 
due to the synchrotron radiation spectrum of a single electron 
(the so-called line of death of synchrotron emission). How- 
ever, at least several cases of GRBs with a < -1/3 have been 
detected (Crider et al. 1997; Preece et al. 1998; Ghirlanda 
et al. 2002, 2003). Second, the typical slope of the low- 
frequency part of the GRB spectrum is a = (Kaneko et al. 
2006), which is not a natural slope of synchrotron radiation 
from shock-accelerated electrons. Third, the high frequency 
spectrum has a highly variable slope /3, not easily explained 
by synchrotron radiation from shock accelerated electrons 
that should produce a fairly standard spectrum with a ^ 1 (as 
in the X-ray afterglows). Finally, for the prompt emission to 
be efficient, electrons are expected to cool fast (Ghisellini et 
al. 2000), but the typical slope of cooling electrons (a = 1 /2) 
is not observed in GRB spectra (Kaneko et al. 2006). 

These well known difficulties of synchrotron emission have 
favored the study of various alternatives, including quasi- 
thermal Comptonization (Ghisellini & Celotti 1999), jitter ra- 
diation (Medvedev & Loeb 1999; Medvedev 2000; Workman 
et al. 2007; Morsony et al. 2008), bulk Compton (Lazzati 
et al. 2000); and photospheric emission (Meszaros & Rees 
2000; Meszai-os et al. 2002; Rees & Meszaros 2005; Pe'er 
et al. 2005, 2006; Giannios 2006; Giannios & Spruit 2007; 
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Pe'er et al. 2007; Thompson et al. 2007; Ryde & Pe'er 2009, 
Lazzati et al. 2009). In this paper we expand the analysis 
of photospheric emission as a GRB prompt radiation mech- 
anism, studying the conditions under which the photosphere 
can produce a high-frequency component with a non-thermal 
power-law shape up to tens of MeV. Photospheric radiation 
has two great advantages with respect to synchrotron emis- 
sion: it does not require a dissipation mechanism and it natu- 
rally provides high efficiency. Photospheric radiation does not 
require a dissipation mechanism if the radiation is released be- 
fore the full acceleration of the fireball and therefore at a stage 
when the fireball still contains a large fraction of internal en- 
ergy. Numerical simulations have shown (Lazzati et al. 2009, 
hereafter LMB09) that this is indeed the case for long duration 
GRBs with a massive star progenitor. The photospheric effi- 
ciency of a typical GRB was found to be in good agreement 
with the observations (LMB09; Zhang et al. 2007). The ob- 
vious weakness of photospheric radiation is that it is custom- 
arily assumed to be thermal, therefore lacking the prominent 
non-thermal tails observed in GRBs. 

The possibility of adding non-thermal tails to the photo- 
spheric spectrum has been investigated before. Pe'er et al. 
(2006) showed that continuous dissipation and/or internal 
shocks at moderate optical depths give rise to a Comptonized 
spectrum with a flat energy spectrum vF{v) oc v^\ Giannios 
(2006, see also Giannios & Spruit 2007) performed a sim- 
ilar study and reached analogous conclusions, focusing on 
magnetic reconnection as the dissipation mechanism. Here 
we argue that a dissipation mechanism that manages to be so 
continuous as to keep the electron temperature at a constant 
equilibrium value under the intense cooling of IC scattering 
is very unlikely. A more realistic assumption is that sub- 
photospheric electron heating takes place as a series of one or 
more episodes of impulsive acceleration and subsequent cool- 
ing. We show that in this case a non-thermal high-frequency 
tail is produced, characterized by a slope that is sensitive to 
a combination of various parameters and which is therefore 
expected to be highly variable during the prompt emission of 
GRBs. 

This paper is organized as follows: in § 2 we discuss the 
origin of the non-thermal tails, in § 3 we discuss the propa- 
gation of the non-thermal spectrum in the optically thick sub- 
photospheric plasma, and in § 4 we present Monte Carlo cal- 
culations of the model. We summarize and discuss our results 
in §5. 

2. SUB-PHOTOSPHERIC DISSIPATION AND SPECTRUM 
FORMATION 

If dissipation takes place in the highly optically thick 
phase of the GRB outflow, multiple scattering and absorption- 
emission processes are likely to thermalize the radiation and 
the baryons, producing a blackbody radiation spectrum at the 
same temperature of the baryons (e.g. Giannios & Spruit 
2007). A more interesting case is the possibility of dissipa- 
tion in the sub-photospheric zone, i.e. for optical depths of 
several. The sub-photospheric zone of a long-duration GRB is 
characterized by the domination of photons over baryons and 
leptons, a condition that, as we will see, is fundamental for 
producing power-law tails with variable slope. In a spherical 
or conical fireball expanding without dissipation (for which 
r = TqR/Rq), the comovingjbaryon density below the satura- 

^ In this and in the next section all quantities are implicitly assumed to be 
in the local comoving frame of the fireball. 



tion radius reads: 
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where Miso is the isotropic-equivalent mass loss rate of the 
central engine, Ri) is the nozzle radius, i.e., the radius at which 
r = 1, and R is the radius of the fireball. The photon density, 
assuming the spectrum is a pure blackbody at thermal equi- 
librium with the baryons, is given by: 
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where Liso is the isotropic-equivalent luminosity of the central 
engine and a = 7.56 x 10"'^ erg cm"^ K""* is the radiation den- 
sity constant. The ratio of photon to baryon densities below 
the saturation radius is therefore independent of distance: 
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where Too. 3 is the asymptotic Lorentz factor of the fireball 
in units of lO"', Rqj is the nozzle radius in units of 10^ cm, 
and Liso.52 is the isotropic-equivalent luminosity in units of 
10^^ erg s"'. Beyond the saturation radius, both the baryons 
and photons densities scale as R~^ (since the photon tempera- 
ture scales as R~^^^, e.g., Meszaros & Rees 2000) and the ra- 
tio is again constant out to the photosphere, where baryon and 
photons decouple. In any reasonable condition, the photon 
to baryon number density ratio of a GRB fireball is therefore 
in the hundreds of thousand. In the absence of pairs, this is 
also the ratio of photons to leptons. If pairs are present, the 
ratio can be altered, but the fact that there are more photons 
than leptons holds since only a fraction of all photons have 
enough energy to create an electron-positron pair (see below 
for more details). Numerical simulations of dissipative fire- 
balls, where the dissipation is provided by the interaction of 
the outflow with the progenitor star, confirm this result. For 
example, LMB09 found a ratio of photons over electrons of 
^ 10^ at the photosphere in their simulation of a typical GRB 
jet with a massive stellar progenitor 

Under such conditions, a non-thermal spectrum can be gen- 
erated by inverse Compton scattering of the thermal photons 
off mildly relativistic electrons. The mechanism of the for- 
mation of the non-thermal tail is the following. Consider a 
dissipation mechanism that produces a population of mildly 
relativistic electrons. For simplicity let us consider a popula- 
tion of thermal electrons with a typical Lorentz factor jg. The 
mechanism can heat repeatedly the same electrons but it is not 
continuous, so that the electrons that are heated are allowed 
to cool before being re-heated. An example of such a mecha- 
nism is repeated shocking by weak shocks, such as those seen 
in LMB09. If the photons dominate over the electrons, the 
electron cooling is very fast and most of the photons never 
scatter off hot electrons, since by the time they scatter off an 
electron it has akeady been cooled through IC interactions 
with other photons. As a consequence, the peak frequency of 
the spectrum is not changed by the dissipation and subsequent 
IC process. However, a fraction (^ < 1 of the photons scatter 
off hot electrons and their energy is increased on average by a 
factor 4/37^. An even smaller fraction of photons C^ ^ 1 is 
scattered twice by hot electrons and its energy is augmented 
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by (4/37^?)^. In this way, a non-thermal tail is built, connected 
to the peak of the original blackbody spectrum. 

To compute the spectral index /3 of the high frequency tai0 
we proceed as follows. Let n^ be the comoving density of 
photons in the primary blackbody spectrum where the popu- 
lation of relativistic electrons is embedded. The peak of the 
F{i^) spectrum is proportional to «^, since the total energy 
density in photons is e ^ n^hvp^ and the peak of the Fiiy) 
spectrum is fpk(i^) ~ e/hfpk = n^. 

The number of photons scattered is proportional to the num- 
ber of relativistic electrons times the number of scattering re- 
quired to cool an electron. Since the energy extracted per scat- 
tering is on averagqj 5E = 4-/3^^hi>, the number of scatterings 
required to cool an electron is 
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Following the above reasoning it's easy to show that 
F{i^i) oc «e«scatt where tie is the electron density and z^i is the 
frequency of the typical blackbody photons that have been 
scattered once vi = A/Jt^^v^^^. The power-law spectral slope 
is therefore given by: 
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Figure [T] shows the value of the spectral index /3 as a func- 
tion of the photon to electron ratio for a sample of combi- 
nations of the remaining two free parameters: the comoving 
peak of the blackbody spectrum /zi'pk and the typical Lorentz 
factor of the accelerated electrons. In all cases the spectral in- 
dex grows logarithmically with the ratio of the photon to elec- 
tron densities. Therefore, outflows very rich in the radiation 
component will have steeper power-law tails. A steepening of 
13 is also observed for increasing comoving peak frequency. 
Finally, more energetic electrons give harder tails, at the price 
of carrying some bumpiness to the spectrum (see below). 

Figure [T] also shows that there is a limit to the photon to 
electron density ratio for which this mechanism is applica- 
ble. As we outlined above, the derivation is based on the as- 
sumption that the probability of a photon scattering off a hot 
electron is less than unity. Therefore, the condition of appli- 
cabihty reads: rij/rie > Wscatt or 
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where e^ and e^ are the energy densities in photons and elec- 
trons, respectively. 

2.1. Repeated dissipation events 

We have so far assumed that there is a single dissipa- 
tion event that energizes the electrons that subsequently cool 

' Note that /3 is defined as the slope of the F(v) spectrum (F(v) oc u^P), 
differently from the Band /3Baiid that is defined as the slope of the photon 
spectrum N(v). Therefore /3 = —(P^^^i + 1). 

* Note that in this analytical derivation we use equations that are formally 
valid in the 7(, — > 00 limit, even though we also consider trans-relativistic 
electrons. In the following, a Monte Carlo code will be used to compute 
spectra more accurately. 



off the blackbody spectrum, producing a non-thermal high- 
frequency tail. It is however possible that the dissipation is 
intermittent, and the spectrum is processed by more than one 
population of hot electrons. For example, in the simulation 
of LMB09 at least 6 shocks are detected between the tj =A 
region and the photosphere. 

If all re-energizations are the same and the electrons always 
have the same typical Lorentz factor 7^, then the process re- 
peats identically and repeated accelerations have the same ef- 
fect as a lower photon to electron ratio. The spectral slope is 
therefore modified as: 
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where «acc is the number of acceleration events. 

If the different energization events produce electrons with 
different typical Lorentz factor -fe, and there is one energiza- 
tion that clearly dominates over the others, then the resulting 
spectrum is expected to be dominated by that energization. 
If, however, the energizations are similar but not identical to 
each other, is is hard to make any quantitative prediction. In 
§ 4 we show with Monte Carlo calculations some examples of 
resulting spectra. 

2.2. Magnetically dominated outflows 

If the outflow is either unmagnetized or moderately magne- 
tized so that the energy density of the magnetic field is lower 
than that of the radiation, the fact that the cooling is dom- 
inated by IC scattering is straightforward. Also in the case 
of Ub > Ui-ad, however, the cooling can be dominated by IC 
interactions if the peak of the synchrotron spectrum is self- 
absorbed (GhiselHni et al. 1998; Giannios 2008). Let eg be 
the ratio between the magnetic energy density and the radia- 
tion energy density eg = Us/Urad- The intensity of the field is 

therefore 

2 

(8) 
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where k is Boltzmann constant. The synchrotron self- 
absorption coefficient for a thermal electron distribution is 
given by a^ = jyC^/ilv^kT) (Rybicki & Lightman 1979), 
where the emission coefficient yV at the synchrotron peak fre- 
quency is y^^^„ = Q.Q'ini.e^B/imgC^). At the synchrotron peak 

frequency i^syn = Q.Ql'y^eB/(mc), the self-absorption coeffi- 
cient reads: 
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where kT = •jgrn^.c^ has been used in the right hand term. 
Combining Eq.|9]with Eq.[8] we obtain: 
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where l'^^^ j^^y is the peak of the thermal spectrum in units of 
keV. The optical depth for synchrotron self-absorption at the 
synchrotron peak frequency Tsyn = a,y^ ^AR therefore reads: 
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where n^AR = Oj at the photosphere. Equation [TT] gives an 
optical depth much larger than unity for moderate values of 
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7e and for eg not much larger than one, showing that our 
mechanism is applicable even to magnetically dominated out- 
flows with eg ~ 100. High frequency power-law tails from 
synchrotron radiation are not expected given the absence of 
non-thermal electrons. 

2.3. Outside the photosphere 

Even though we have derived the spectrum and applicabil- 
ity conditions for sub-photospheric dissipation, the unique re- 
quirement for this mechanism to work is that the electrons are 
cooled by IC scattering with only a fraction of the photons. 
This can naturally take place in the optically thin part of the 
spectrum. Outside the photosphere the number of IC scatter- 
ings per photon is < 1, however the number of IC scatterings 
per electron is still large, the condition for applicability of this 
derivation outside the photosphere reads: 



Tt > 




(12) 



which, according to Eq. |3] can be a substantially small num- 
ber. Therefore, the non thermal spectrum can be generated 
with this mechanism out to large radii, approximately 100 
times bigger than the photosphere. This has important con- 
sequences for the extent of the spectrum at high frequencies. 
If the episodic dissipation that we envisage here can take place 
at large radii, well outside the photosphere, the radiation up to 
GeV frequencies could be accounted for 

3. RADIATION TRANSFER IN THE COOLED ELECTRONS 

The high frequency power-law tail is produced in a region 
of Tt ~few. In a relativistically expanding medium, all pho- 
tons undergo a number ^ tt of scattering off cold electrons 
before leaving the region (e.g., Pe'er et al. 2005). Since the 
change of energy of a photon scattering off a cold electron is 
on average 61^ /i^ ~ hv/nieC^, a cutoff in the spectrum at a co- 
moving photon energy hv = 51 I/tj- keV is created. As long 
as the Thomson depth is not too high, the power-law tail is 
preserved up to several tens of keV, likely extending into the 
tens of MeV in the observer frame. 

Another effect of photon propagation, once the non-thermal 
tail has been added to the blackbody spectrum, is the creation 
of electron-positron pairs via photon-photon collisions. If the 
pairs created in such a way do not outnumber the original elec- 
trons, the only effect of their creation is the absorption of all 
the photons above hv = 5l\ keV. This would introduce a cut- 
off ?A hv = 511 keV which is irrelevant since the downscat- 
tering off cold electrons discussed above introduces an even 
more severe cutoff in the spectrum. If, instead, the electron- 
positron pairs outnumber the original electrons, a pair pho- 
tosphere is created at a larger radius (Pe'er & Waxman 2004; 
Rees & Meszaros 2005). Typically, the pair photosphere is lo- 
cated a factor 3-10 farther than the electron photosphere. The 
processes described in section 2 do not change if the opacity 
is provided by pairs rather than by electrons. The number of 
pairs, however, is not independent from the number of pho- 
tons, since the pairs are created by photon collisions. The 
photon to pair density ratio can therefore be easily calculated: 
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where f3' is the slope of the spectrum that created the pairs 
and needs not be equal to /?, the slope of the spectrum created 



once the pair photosphere has been established. The spectral 
slope f3 now reads: 



Log 



Pe+e ■ 



4 7," f hi^pk \ 
3 7^.-1 \ nigC^ J 



1-/3' 



Log(4/372) 



(14) 



If, however, the cyclical process of creating the power-law 
tail, producing pairs from the photons above threshold, pro- 
ducing a new tail etc. takes place several or more times, it 
is likely that the spectrum will saturate such that f3e+e- ^ f3', 
resulting in a saturated slope: 
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Figure|2]shows the value of the spectral slope as a function of 
the peak frequency for a range of the electrons' Lorentz factor. 
Interestingly, in the vast majority of cases an episodically dis- 
sipative pair photosphere yields a spectral slope 1 < /3sat < 1 .5 
in fairly good agreement with the observational results of 
BATSE (Kaneko et al. 2006). 

4. MONTE CARLO SIMULATIONS 

In the past two sections we have discussed the radiation 
properties in the comoving frame of the GRB outflow. In 
this section we move our reference frame to the observer and 
present Monte Carlo (MC) simulations of the radiation in the 
observer frame. In this frame, the spectrum depends on an ad- 
ditional free parameter: the bulk Lorentz factor of the fireball 
r. The MC code assumes that a blackbody spectrum of pho- 
tons is embedded in a relativistically expanding slab of elec- 
trons with a certain Thomson opacity rj-. At the beginning 
of the simulation, the electrons are accelerated to a unique 
value of their random Lorentz factor 7^ . The position of the 
photons, ten million per simulation, is initialized as a random 
location in the slab, and after each scattering it is updated ac- 
cording to the relevant distribution of traveled distances: 
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where / is the distance traveled by a photon between scatter- 
ings and /() = Tr ' is the mean free path of low frequency pho- 
tons. The full Klein-Nishina cross section was used to com- 
pute lo as a function of the photon energy. At each scattering, 
the energy of the photon is modified in two possible ways. If 
the electron off which the scattering takes place is relativis- 
tic {{je- 1) > hv/ifigC^), the new photon energy is generated 
according to the IC energy spectrum for a single scattering re- 
ported, e.g., in Rybicki & Lightman (1979, Eq. 7.24). The 
Lorentz factor of the electron is also updated to take into ac- 
count the energy given to the photon. If, on the other hand, the 
electron is non-relativistic, Compton scattering is assumed, 
and the new photon energy is computed from the Compton 
equation for a randomly generated scattering angle (e.g., Ry- 
bicki & Lightman 1979, Eq. 7.2). Also in this case the energy 
of the electron is updated accordingly to the energy loss of the 
photon. 

Adiabatic energy losses are not included in the code be- 
cause the cooling time scale of the electrons is much faster 
than the adiabatic time scale. The comoving cooling time 
scale of electrons is 
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while the comoving adiabatic time scale 
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phot 
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is of the order of seconds. Another effect that is not consid- 
ered is the cooling of the photons as a consequence of the 
expansion of the electrons (adiabatic cooling of the photons). 
Due to the small optical depths considered, this effect is also 
negligible. 

At each step the position of the photons is checked and all 
the photons that have reached the surface of the slab are col- 
lected in the output spectrum. The code is stopped when 1/3 
of the photons are collected outside of the slab. The scattering 
is therefore computed well beyond the step at which all the 
electrons have cooled in order to properly take into account 
the radiation transfer into the cold electrons. The output spec- 
trum is then processed to take into account the fact that the 
emission is produced on a spherically expanding surface. 

Figure [3] shows result of three MC runs for a GRB with 
observed peak frequency /ii'obs ^ 500 keV and Lorentz fac- 
tor r = 1000. The three runs are performed for dissipation at 
tt = 2 and with different values of the photon to electron den- 
sity ratio and typical electron Lorentz factor 7^,. The figure 
shows that indeed a prominent power-law tail is added to the 
thermal spectrum. The tail extends to at least several tens of 
MeV in the observer frame, in good agreement with observa- 
tions. The main difference with the predictions of § 2 is that 
the power-law slope is steeper by approximately 0.5. The pre- 
dicted slopes are /3 = 1.23 (solid line), /? = 1.53 (dashed line), 
and /3 = 1.98 (dotted line). However, the slopes measured in 
the figure are /3 = 1.72, /3 = 2, and /? = 2.5, respectively. The 
reason for the discrepancy is that the equations of § 2 assume 
that the electrons remain at 7 = 7^ until they cool and disap- 
pear. In reality, photons that scatter more than once off a hot 
electron are likely to find the electron slightly cooled at the 
second scattering. This has the effect of producing a slightly 
steeper power-law slope. As we will see in the following, this 
can be cured by assuming that there are more soft photons (so 
that not all scattering extract the same energy from the hot 
electrons). Figure |4] shows instead a comparison between the 
solid line of Fig.[3]and two analogous simulations in which the 
dissipation is assumed to take place at a higher optical depth 
Tt = S and tj = 16. As expected, the three spectra are virtu- 
ally identical at low frequencies, but the rj = 8 spectrum has 
a cutoff at ^ 40 MeV and the tt = 16 spectrum has a cutoff at 
^ 10 MeV due to direct Compton scattering. 

Figure |5] shows a comparison of spectra with analogous 
characteristics but with different observed peak frequencies 
and bulk Lorentz factors. The figure confirms that the mecha- 
nism can reproduce spectra with both low and high peak fre- 
quency. An interesting feature that can be seen in Fig.|5]is the 
bumpiness of the dotted spectrum. It turns out the spectrum 
looks bumpy for jg > 3, i.e., when the photon energy gain per 
scattering is larger than an order of magnitude. 

4.1. Repeated dissipation events 

If the spectrum is produced by multiple dissipation events 
with different characteristics (different value of 7^), the ana- 
lytic predictions do not work and only a numerical calculation 
can give the appearance of the final spectrum. Figure|6] shows 
the outcome of two possible situations, one in which there 
is a dissipation event which dominates and another in which 
the three dissipation events are comparable in strength (and, 
therefore, ■^e)- 



The solid line shows the results of a MC simulation in 
which three energization events are considered: a small one, 
yielding 7f=2; an intermediate one, yielding 7(. = 5 and a fi- 
nal small one, yielding 7^ = 2. As qualitatively predicted, in 
that case the spectrum is dominated by the strong event. The 
main effect of the smaller events is to pre-broaden the input 
photon population, allowing for a smoother power-law rather 
than a bumpy one (compare to the dashed line in the same 
Figure |6l). The same conclusion holds even if the three en- 
ergizing events are characterized by fairly similar energetics: 
the resulting spectrum matches fairly well the prediction for a 
single energization with the largest 7;,. 

4.2. Non-thermal input spectrum 

We have so far considered the spectrum arising from the 
Comptonization of a primary thermal spectrum. However, ob- 
served GRBs not only have non-thermal high frequency tails, 
they also have non-thermal low-frequency tails, with a typi- 
cal a = (where, again, F{v) oc j^~"; Kaneko et al. 2006). 
Non-thermal low-frequency power-law tails require a decou- 
pling of the photons from the electrons below the photosphere 
and studying the origin is beyond the scope of this paper In 
any case, it is interesting to see what is the consequence of 
a non-thermal low-frequency spectrum on the Comptonized 
high-frequency tail. 

To investigate this issue, we performed MC simulations 
with thermal spectra modified to have a v^ and i/' low- 
frequency behavior We find that broadening the input spec- 
trum has two important consequences. First, the slope of the 
Comptonized spectrum obtained is in better agreement with 
the prediction of Eq. |71 As seen in Fig. |7] the i^' input spec- 
trum already hardens the power-law tail, while the j/" input 
spectrum provides an very good agreement with the predic- 
tion (shown with a thin dashed line). Figure [8] shows that be- 
sides the better agreement with the theoretical slope, allowing 
for a flat low-frequency input spectrum removes the marked 
oscillations (or humps) seen in the high-frequency tails in sit- 
uations where the energy gain per scattering is large, again 
improving the agreement between the spectra in the figure and 
the observations. 

4.3. Non-thermal relativistic electrons 

We have so far assumed that the energy dissipation pro- 
duces a thermal population of relativistic electrons character- 
ized by a unique Lorentz factor 7^,. In some cases, however, 
dissipation can lead to the acceleration of non-thermal elec- 
trons, with typical energy spectra dng/d^ oc 7"'', with p >2. 
Whether or not the presence of the non-thermal electrons af- 
fect the spectra discussed above depends on the slope of the 
electron distribution and on the fraction of electrons that are 
contained in the non-thermal tail. The strongest influence 
on the spectrum is observed if all the electrons are acceler- 
ated in the non-thermal tail with a hard typical slope (e.g., 
p ^ 2.3). In that case, the photons that IC scatter off hot 
electrons produce a steep non-thermal spectrum with slope 
13 = (p- l)/2 - 0.65, 13 = p/2= 1.15 OT (3 = 0.5, depending 
on the cooling regime. Since the majority of the electrons 
do not scatter off hot electrons but rather off already cooled 
electrons, the power-law tail does not connect smoothly with 
the primary spectrum. The results of MC simulations of this 
scenario are shown in Figure |9] If the primary spectrum is 
assumed to be a pure blackbody, the Comptonized one has 
a prominent power-law tail that does not connect smoothly to 



Lazzati & Begelman 



the primary photons, creating a hardening of the spectrum that 
is not usually observed in GRBs (e.g. Kaneko et al. 2006). 
Only if a very soft input photon spectrum is assumed is a 
Band-like spectrum obtained (Fig. |9]l. It is clear, however, 
that the mechanism described in this paper produces spectra 
in better agreement with observations if the electron accelera- 
tion events contain only a moderate amount of energy and do 
not produce a sizable non-thermal tail in the electron popula- 
tion. 

5. DISCUSSION 

We have presented a method by which a primary blackbody 
spectrum of a GRB outflow photosphere can be Comptonized 
into a non-thermal high-frequency power-law tail. In this sce- 
nario, the bulk of the GRB prompt emission is thermal pho- 
tospheric and the non-thermal tail contains a relatively small 
fraction of the radiation energy. The majority of the energy 
of the radiation is therefore released before it is converted 
into bulk kinetic energy, solving the problem of identifying 
a dissipation mechanism with high efficiency (LMB09). The 
non-thermal tail is produced by thermal electrons accelerated 
to mildly relativistic energies by an intermittent dissipation 
mechanism. The electrons are allowed to cool through IC in- 
teractions with the photon field between acceleration events. 
This dissipation mechanism is not required to have an effi- 
ciency larger than a few per cent, since the power-law tail 
only contains a fraction of the total energy of the radiation. The 
spectral slope of the high energy radiation depends on the typ- 
ical Lorentz factor of the accelerated electrons, on the peak 
frequency of the thermal photon spectrum, and on the ratio of 
the photon to lepton density. Differently from previous work 
(Pe'er et al. 2005, 2006; Giannios 2006) we assume that the 
electrons are accelerated intermittently and not continuously 
and that the electrons are allowed to cool among acceleration 
events. In addition, we assume that the energy density in radi- 
ation is larger than that in the relativistic electrons, a situation 
that is easily realized in jets born inside massive progenitor 
stars (LMB09). 

Compared to the standard optically thin synchrotron, this 
mechanism can explain the steep low-frequency slopes ob- 
served in the early phases of some bright bursts (Crider et al. 
1997; Preece et al. 1998; Ghirlanda et al. 2002, 2003) and 
the transient thermal bumps detected is several events (Ryde 
2005; Ryde & Pe'er 2009; Ryde et al. 2010). Bursts for 
which no dissipation takes place should not display any non- 
thermal features. A complete discussion of the low-frequency 
non-thermal tail is however beyond the scope of this paper 
and will be presented elsewhere. Explaining a typical low- 
frequency tail F{i/) oc v'^ is indeed a challenge for any model 
and not only for the photospheric scenarios presented here 
and elsewhere (Meszaros & Rees 2000; Meszaros et al. 2002; 
Rees& Meszaros 2005; Pe'er et al. 2005, 2006; Giannios & 
Spruit 2007; Pe'er et al. 2007; Thompson et al. 2007; Ryde 
& Pe'er 2009, LMB09). 

Radiation from dissipative photospheres has been investi- 
gated before (Pe'er et al. 2005, 2006; Giannios 2006; Gi- 
annios & Spruit 2007; Beloborodov 2010). Giannios (2006) 
and Giannios & Spriut (2007) concentrate on the case of a 
continuous energy injection or "slow heating", different from 
our impulsive acceleration assumption. Pe'er et al. (2005, 
2006) consider instead two possible scenarios: "slow heat- 
ing" and impulsive heating by internal shocks. They find that 
the bulk of the photons are shifted in energy once a steady 
state electron population is attained. They calculate that, for 



a fairly wide parameter range, the comoving peak frequency 
is of few tens keV. Attaining such a steady state population 
requires however a large optical depth and as a consequence 
we do not find the same results. To understand the reason for 
this, consider an injection episode at an optical depth of a few. 
Due to IC scattering, the hot electrons are efficiently cooled 
and their energy given to a small fraction of photons that are 
shifted to high frequencies. For small optical depth, the en- 
ergized photons reach the photosphere with only negligible 
losses and decouple from the electron population. The en- 
ergy given to the electrons by the dissipation event is therefore 
given to a small number of photons that produce a power-law 
tail as shown in the figures. In order to obtain a steady state 
electron population able to exchange energy effectively with 
the radiation field, a large number of scattering per photon is 
required, so that the high-frequency photons can return their 
energy to the electron population. Large optical depths are 
therefore required for electron and photon populations to set- 
tle in a steady state configuration dictated by the balance of 
heating and cooling. Such optical depths are much larger than 
the ones envisioned in this work. In any case, a fraction of the 
dissipated energy is left in the electron population that ther- 
malizes at a temperature higher than the one before the dissi- 
pation. For the parameters adopted in this paper the effect of 
the increased electron temperature is negligible because the 
increase of the electrons temperature is very small for the op- 
tical depths considered. Finally, Beloborodov (2010) inves- 
tigated the radiation produced by a well-defined dissipation 
mechanism: nuclear and Coulomb collisions within a bary- 
onic outflow with a substantial population of neutrons (see 
also Beloborodov 2003; Rossi et al. 2006). He finds that a 
typical spectrum with a high-frequency slope Fiiy) oc v~^'^ is 
obtained. Differently to our result, the non-thermal spectrum 
is due to the presence of non-thermal leptons. 

A characteristic of any photospheric model is that the radi- 
ation is released at a small radius rph ^ 10'""'-' cm (e.g. Rees 
& Meszaros 2005). The main consequence of such a small 
radius is that the compactness of the region is large (Pe'er 
& Waxman 2004) and therefore photon-photon interactions 
result in the production of electron-positron pairs. This can 
have two consequences. If the newly created pairs do not 
outnumber the original electrons, the only consequence is the 
presence of a sharp cutoff in the spectrum at hv = 511 keV 
in the comoving frame (^ 100 MeV observed). If, on the 
other hand, the pairs outnumber the original electrons, a new 
photosphere is created, typically at a radius 3 to 10 times 
larger than the original photosphere (Rees & Meszaros 2005). 
The scenario we developed in this paper is not dependent on 
whether the photosphere is due to the original electrons or 
to pairs, the only difference being that a pair photosphere 
has a lower photon to lepton ratio and produces therefore a 
steeper high-frequency spectrum. In addition, since the num- 
ber of pairs and the number of photons are not independent, 
we showed that a pair photosphere with at least several ener- 
gization events would produce spectra with a typical slope 
1 < /3sat < 1-5, in good agreement with observations (e.g., 
Kaneko et al. 2006). 

The Fermi satellite has recently shown that at least some 
bright GRBs have spectra extending well into the Ge V regime 
(Abdo et al. 2009abc; 2010; Ghisellini & Ghirlanda 2010; 
Granot et al. 2010). For any reasonable combination of pa- 
rameters, GeV radiation cannot be produced at the photo- 
sphere, since it would be immediately absorbed by photon- 
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photon colHsions to produce pairs. In this scenario, the GeV 
emission can be either internal, i.e. produced by dissipation 
above the photosphere (see § 2.2; Toma et al. 2010) or ex- 
ternal, due to the interaction of the outflow with the interstel- 
lar medium (Kumar & Barniol Duran 2010; Ghirlanda et al. 
2010; GhiselUni et al. 2010). 

In terms of predictive power, the most notable prediction 
of this model is the presence of substantial variability of the 
high-frequency spectrum due to the dependence of the spec- 
trum on parameters that are expected to be highly variable 
such as, e.g., the peak frequency of the thermal spectrum 
in the comoving frame (LMB09). Such variability might 
be however hard to detect since it is expected to have short 
time-scales and would therefore be averaged out in a time- 
integrated spectrum. The dependence of the spectral slope on 



the peak frequency of the spectrum opens also the possibility 
of correlations between the peak frequency and the spectral 
slope at high frequencies. Since, however, the slope depends 
on two more parameters (the ratio of photons to electrons and 
the typical Lorentz factors of the electrons) it is not obvious 
that such correlations should be visible in experimental data. 
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Fig. 1. — The high energy spectral index /3, defined through F{y) onv ^ plotted versus the ratio of photon to electron (or electron-positron pairs) density. 
Several combinations of the parameters 7^ and Az/pk are shown. The dark band shows the range of ratios observed at the photosphere of the LMB()9 simulation. 
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Fig. 2. — Saturated slope of the high frequency tail of the spectrum for a pair photosphere after a few cycles of power-law spectrum generation and pair 
production. 
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dashed line shows the predicted slope for the strongest dissipation (7 = 5; /3 = 1.06). 



14 



Lazzati & Begelman 



10' 



T3 
0) 

_N 

Iq 10"^ 



o 

c 

S 10"' 






10 




r = 1000; /ii^„fe,=500 keV; r^^ =2 
7, =2; n /n, =3 xlO^ 



Black body 
.1 



10^ 

/iz/(keV) 



Fig. 7. — Comparison of Monte Carlo spectra for different primary spectra. Tlie solid line is the same spectrum shown in Figure |3] The figure shows how 
broadening the low-frequency part of the input spectrum results in a better agreement with the slope predicted theoretically in § 2. 



Non-thermal GRB photospheres I 



15 




10^ 

hu (key) 

Fig. 8. — Comparison of Monte Carlo spectra for different primary spectra. Tlie solid line is the same spectrum sliown witli a dotted line in Figure |5] The 
figure shows how broadening the low-frequency part of the input spectrum results in a smoother power-law, without the humps caused by the large energy gain 
per scattering. 



16 



Lazzati & Begelman 



10° 




1 1 1 1 










* ' • (^^^ ^^. 






-•^'''^ ' * ^\. 


■ 




^^^ * * ' ^x 






y/^ '♦ ^\. 




^"^^ 


y^ ' ♦ ''• ^ ^v 




-o 


./^ * ' ^ \. 




(U 


.^ ' " ■ ^ ^v 




_N 


/^ ' ' • ''. ^ ^\. 




Id 10"^ 


/ ' ' • •''-^\ 


- 


E 


/ "*^^ 




o 




\*^"»^ 


c 




\ "- < .. 


^ 






ro 








r = 1500; hiy^^^ =1000 keV; Tt=2 
7, =3; n^/n,= lO^p =2.3 


\ ^ 

\ ^ 

\ s 


^ 






^ 




\ 




Thermal electrons 


\ 




Black body 


>y^ 


10-^ 


- u' u' 


\ 




1 ■ ■ 1 


\ ^ 



10' 



10^ 

huikeV) 



10^ 



Fig. 9. — Monte Carlo spectra with non-thermal electron population. 



